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Abstract 

A dense Bose gas with hard-core interaction is considered in an optical lattice. We study the 
phase diagram in terms of a special mean-field theory that describes a Bose-Einstein condensate and 
a Mott insulator with a single particle per lattice site for zero as well as for non-zero temperatures. 
We calculate the densities, the excitation spectrum and the static structure factor for each of these 
phases. 
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I. INTRODUCTION 



Optical lattices have opened an exciting field of physics. We are expecting new phenomena 
in comparison with continuous systems due to the lack of Galilean invariance. Among the 
most obvious consequences of the lattice structure is the formation of lattice-commensurate 
ground states like the Mott insulating phase. The latter was observed in experiments 1, 2]. 

Light scattering on a Bose gas is strongly affected by the nature of the quasiparticles |3|,|4j. 
This provides a useful experimental tool to distinguish between different phases of the Bose 
gas. We expect, for instance, that light scattering by the gapless quasiparticle spectrum of the 
Bose-Einstein condensate is quite different from light scattering in the gapful quasiparticle 
spectrum of the Mott insulator (MI) . A physical quantity that is directly related to light 
scattering is the structure factor {4], 3, 0, Q, Is, 9] which was measured in the case of a BEC 
[2I. This is also of interest in the Mott insulating phase 10, 11]. In this paper we will study 
the static structure factor. Of particular interest is its behavior near the transition from the 
BEC to the MI. We will also study the phases and phase transitions for zero temperature 
as well as for nonzero temperature. For this purpose we will consider here a hard-core Bose 
gas in an optical lattice. Our model can be understood as a projection of the more general 
Bose-Hubbard model in the vicinity of those points of the phase diagram, where two adjacent 
Mott lobes meet (Fig. [p. This is similar to the picture which was applied to the tips of 
the Mott lobes in a recent paper by Huber et al. |8j. It is based on the following idea. The 
number of bosons per site is fixed in the Mott state. For adjacent Mott lobes this means 
that the corresponding Mott states differ exactly by one boson per site. Now we consider 
two adjacent lobes with n and n + 1 (n > bosons per site), respectively and assume that 
the chemical potential is fixed such that the ground state is the Mott state with n particles 
per site. Low-energy excitations of this state for a grand-canonical system are states, where 
one or a few sites (let's say k > 1 sites) have n+1 bosons, all other sites have n bosons. The 
k excessive bosons are relatively free to move from site to site on top of the n Mott state. 
Therefore, the physics of these excitations can be described approximately by the tunneling 
of the k excessive bosons alone. Due to the repulsion of order U, assumed to be not too 
small, it is unlikely that a site with n + 2 bosons is created. Consequently, these excessive 
bosons form a hard-core Bose gas. 

A hard-core boson can be represented by a pair of locally coupled spin-1/2 fermions. Here 
we will use the model introduced to study the dissociation of bosonic molecules into pairs of 
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spin-1/2 fermionic atoms 12j. A functional integral with a two- component complex field was 
derived for this model that allowed a mean-field approximation. The latter revealed a zero- 
temperature phase diagram for a grand-canonical enesemble of bosons with three phases: an 
empty phase, a BEC and a MI. This result is remarkable, since previous mean-field calcula- 
tions, using the same type of fermionic model as starting point, (e.g. using an N — > oo limit 
[13I ) did not give the entire phase diagram by ignoring the MI phase. Aparently, the form 
of the phase diagram depends crucially on the type of Hubbard-Stratonovich transformation 
that replaces the fermionic (i.e. Grassmann) fields by bosonic (i.e. complex) fields. 

The structure of this paper is as follows. In Sec. 2 we discuss briefly the model and some 
related physical quantities. A mean-field approximation is used in Sec. 3, where we derive the 
phase diagram. Gaussian fluctuations around the mean-field solution and its consequences 
are considered in Sec. 4. In this section we calculate the static structure factor. And in Sect. 
5 we discuss our results. Details of our calculations are given in the Appendices A and B. 

II. THE MODEL 

A hard-core Bose gas is considered on a d- dimensional hyper-cubic lattice with M sites. 
Bosons can tunnel between these sites from a lattice site r to any of the nearest-neighbor 
sites r' . Each boson consists of two fermions. Single fermions in this model can not exist 
(i.e., in contrast to the model studied in Ref. 12J we do not allow dissociation of the bosonic 



molecules). The Hamiltonian of the <i-dimesional system is 

H = Yl C It C r't C I| C r'l -/^EE 4a C ra, (1) 

(r,r'> r cr=U 

where the sum in the first term on the right-hand side is restricted to nearest neighbors, and 
c\. a and c ra are the creation and annihilation operators of the fermion with spin a at site r, 
respectively. The first term describes tunnelling of local fermion pairs in the optical lattice 
between nearest-neighbor sites with rate J. The chemical potential fi controls the number 
of particles in a grand-canonical ensemble. The latter is given by the partition function 

Z = Tre- 06 . (2) 



The partition function can also be written in terms of a Grassmann integral 14j as 



- A ^D[M] (3) 



with imaginary time-dependent conjugated Grassmann fields ip, ip and the action 



A I dt 
'o 



J 



(r,r') 



(4) 



Applying the concept of linear response to the quasiparticles, we obtain the dynamic 
structure factor with the definition Q 



S(q, U) = Y I ~ ft) \ m )\ 2 ^ - ^n,m) 



(5) 



n,m 



where is the Fourier transform of the atomic density operator at wave vector q and io n , m 
is the frequency difference between energy level m and n. We treat an excited state n as one 
of the quasiparticles energy level with fojj n $ = e n within our approximation. Then the static 
structure factor is defined as 



%) = f S(q,Uj)duJ = ^{(PqP-q) - \ (pq)\ 2 ). 



(6) 



For small but nonzero T, when quasiparticles can be treated as noninteracting within the 



accuracy of the Bogo 
optical lattice we get 

where 
For T = 



iubov approach for the weakly interacting Bose gas modified for an 

a 

S ( q ) = ^l cot h^, (7) 



1 

9q = 1 - -lY C ° Sqi 



S(q) 



i=l 



Jg q 



(9) 



where e q is the quasiparticle spectrum. 

This result was obtained within the mean-field approximation to the weakly interacting 
Bose gas. It is interesting to compare this result with the result obtained by a direct calcula- 
tion for our model within our mean- field approach. Using Eq.(jn]) the static structure factor 
can be calculated as [3] 

r,r' 

where the truncated density-density correlation function is 



C rr i = (n r n r i) — (n r )(n r i). 



'IF 
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Introducing spatial coordinates p = {r, t} the static structure factor reads 

S W = TTKfYl l dt [ dt'Crtyve^-^. (12) 

P iV r ^ Jo Jo 

This expression will be studied within mean-field theory and Gaussian fluctuations. 



III. MEAN-FIELD APPROXIMATION: PHASE DIAGRAM 

The mean-field approximation cannot be directly applied to the Grassmann fields. There- 
fore, we perform a Hubbard- Stratonovich transformation in order to replace the Grassmann 
fields by complex fields. This leads to an effective action depending on complex fields. As 
we already mentioned in the Introduction, the form of the Hubbard-Stratonovich transfor- 



mation is not unique. Here we use the version given in Ref. 12| because it provides the 
proper phase diagram. First of all we decouple the fourth order term in the action using the 
identity 

eXP I Yd ^ ^\4Wr' I = / [d<f>][dx] exp - §rK)4r> ~ ^ XrXr- 
\ {r,r') I ' L r,r' r 

- y^(i0 r + Xr)VV,lVV,2 - ^(#r + Xr)A,lA,2 



(13) 



Summation over t in the last formula is implicitly assumed. In this formula the expression 

V r y = J ( — 5| r _ r/|> i + 25 r y ) (14) 



2d 

was introduced. Moreover, the introduction of the auxiliary field x makes the matrix v 
positive definite. 

A subsequent integration over Grassmann fields leads to 



^efT 

with 



jT dtj^^u-^ + ^^XrXr-logdetG x | (15) 



, -id> — x d T + a \ A l/2d, r, r' — nearest neighbors 

G = ( _ _ , v = J(cu + 21), Lu r y={ 

d T — fi id) + x I |0, otherwise 

The partition function now reads 

Z = J e- A ^ x) D[ct),x}- (16) 



We can perform a saddle-point integration to calculate physical quantities. From a phys- 
ical point of view we have to minimize our action to get the classical trajectory (i.e. the 
macroscopic wavefunction) of our system. Fluctuations around this trajectory are caused by 
thermal and quantum effects. The mean-field solution characterizes the condensed phase, 
in which \<p\ has a nonzero value. Fluctuations describe quasiparticle excitations above the 
condensate. In order to proceed within a mean-field approximation we assume that quantum 
fluctuations are small. 

Minimization of the action gives us two coupled linear equations between the complex 
fields 4> and x- 

m- = o = , (17) 

X = -i2JG{<P-i X ) 



where 



A solution of these equations is 



C - 1 y 1 

fiZs 10)2/9 + ^ /A + U} 2 n 



1 2ic 
G = -, x 



J 3 

In these equations u n = (2n + 1)tt//3 are the Matsubara frequencies of the fermions. This 
leads to (we perform the rescaling 4|0| 2 /9 — > \<f)\ 2 ) 



J ^ (^ m ^ l ) («) 

and within the saddle-point integration we obtain the expression for the local density 

1 dlogZ = I 1 jl / e /V^V2_i \ 

It should be noticed that this expression for J = and T = gives correct result, which 
can alone be obtained by direct calculations. The corresponding phase diagram is depicted 
in Figs. 1 and 2 with the phase boundary given by <p = 0. 

With the result given in Eq. (l20~l) we can evaluate the local compressibility: 

(21) 

For T = and small |0| 2 it reads 

< 22 > 
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The local compressibility is not divergent wit 



typical in the case of optical lattices [16|, 0, \l8 \. 

We see that this model gives three phases. It can describe both the dilute regime and the 
dense regime, as a consequence we have the empty phase along with the MI phase. 

A few words about the type of the phase transitions. Near the phase transition the order 
parameter is small, so that we may apply a perturbative expansion : 

1 |0| 2 1 1 |0| 2 1 1 '^' 4 



lin our mean-field approximation, which is 



A 



cff 



J 9 



l Y 

R ^ 



(3^ ^ + 9 2/3^(^ + ^)2 81 



+0 



(23) 



can be negative 



>o 



The first term can be negative, such that the order parameter may differ from zero. The order 
parameter vanishes continueously indicating that we have a second order phase transition. 



IV. GAUSSIAN FLUCTUATIONS AROUND MEAN-FIELD SOLUTION 



The complex fields <p x and Xx are expected to fluctuate about the SP solution due to 
thermal and quantum effects. Denote A = icf) + x and A = i<p + x, then 



G- 1 = G - 1 + 



-5A 
5A 



where 



Go 1 



-A d T + fj, 

dr-fJ, A 

Applying the Taylor expansion ln(l + x) = x — x 2 /2 + ... we get 



log det GT 1 = tr In G" 1 = tr In 



G 1 + 1 


( -6A 










K 


8 A J 



« trln G n 1 — -tr 



G r 



-SA 
5A 



Calculating the trace in p = {q, u} representation we get 

Z~ J D[5<j>]exp [-5A eS } 



(24) 



(25) 



(26) 



where 5A e g is given in Appendix A. This has the form of an inverse Green's function G 1 with 
imaginary time for quasiparticles, which is described by fields 8(f). By applying the spectral 



representation of the Green's function we can identify the poles of function G, namely iuj 
with the excitation spectrum of the quasiparticles. The exact solution in the condensed 
phase is 



e q = ^g q (J*-^)+gy. (27) 
In the dilute regime, where fi + J is small, and for small q we get the Bogolubov spectrum 



with an effective mass of the bosons ~ 1/J and sound velocity ~ \f2Jjj, (we used ji = J + 



This has the same form as the quasiparticle spectrum obtained in the 



In the empty phase and in the MI phase we have excitations with a gap: 

e q = - J + Jg q . (28) 

The symmetric form of the results for the MI phase and for the empty phase is due to the 
particle-hole symmetry of our model. 

In the dilute regime (small ft = J + /i) and for low temperatures, including fluctuations, 
we obtain 

E_ f dq , .3 f d 3 q 1 
^ nq ~^Jjir { b ) ~ 2 J (29) 

This result is in good agreement with a weakly interacting Bose gas with an effective chemical 



potential \x [191 ]. 



A. Quantum Fluctuations versus Thermal Fluctuations 

By defining the condensed density in the following way 

n Q = lim (cLcr^c^Cr^) = lim (^^'V'rV'r') 
\i — • r'|— >oo |r— r'\— >oo 

we arrive at the expression for the condensed density, which can be calculated as 
lim / dt / dt' 1 



\r—r'\—>oo 



o Jo \ (-9! + /x 2 /4 - A rt A rt )(-d' t 2 + fi 2 /A - A r , t ,A r , t , 
Then fluctuating fields A = A + 5A, expanding the above expression up to the second order 
of 5 A and using the Green's function given in Appendix A, we arrive at 

n Q = iJ2 ~/ 2) + Sn , (30) 
where the correction due to the quantum fluctuations to the mean-field result is 
6no = ( J2 ~ MV f _d^$k9k + {J 2 - /i 2 ) f J^L Bktk (J 2 - Z^ 2 ) 2 f ddk B l 



J 3 J (2vr) d e k 4J 3 J (2vr) d 4J 3 J (2vr) d t k 
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3(J 2 - /x 2 ) V /• Blg k 
4J 5 7 (27r) d e fc 

The main correction due to the thermal fluctuations are already included in our mean-field 
theory, where the condensed density is given by 

n.= !^, (3D 

and 1 0| 2 can be determined from Eq. (TT9l) . 

The effect of quantum fluctuations and thermal fluctuations is depicted in Fig. HI We see 
that both of them lead to a depletion of the condensate, but the quantum depletion alone 
does not change the transition points. 



B. Static Structure Factor 



The static structure factor is a Fourier transform of the truncated density-density corre- 
lation function. The latter contains nonlocal term 

rP fP 

(^V^Vr^r'^rWr') = dt dt' 



o Jo \ K -A rt A rt - d 2 t + // 2 /4)(-A r / t ,A r , t , - dj + /i 2 /4) 

(32) 

n mean-field approximation the truncated density- density correlation function should vanish 
12j. Taking into account fluctuations, we write again A = Ao + 5 A, A = Ao + 5 A and 
substitute it into Eq. fl32l) . Then expanding up to the second order and Fourier transforming 
the obtained expression, after a direct but lengthy calculations we get for small wave vectors 
q and for low temperatures T in the BEC phase 

Qf \ ~ /i 2 ) Jg q ,,(3e q 

%)~— j^— coth— , (33) 

where n is a total density of particles. 

In the dilute regime, i.e. close to the empty phase, when n ~ (J + fi)/J and J — \i ~ 2 J 
we obtain 

%)~^coth^, 

€q n 

which is in agreement with the well-known result for the weakly interacting Bose gas |5|]. 
In the dense regime, i.e. close to the Mott phase when n ~ 1, the static structure factor 
vanishes. The dependence of the static structure factor on q in the dilute regime is depicted 
in Fig. [5j 
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The static structure factor measures the density- density correlations in q space. Assuming 
that the thermal energy is much larger than the gap (i.e. ksT 3> A), the maximum of the 
structure factor appears at q values for which the following condition holds 

Jg q « A, (34) 

where A = /i — J is an energy gap in the excitation spectrum. Thus fluctuations begin to 



feel each other at the distance equals to their effective size r~l/g~l/ y/fj, — J. 

Spatial correlations can be calculated by Fourier transforming the static structure factor. 
Near the phase transition for T = we have in the condensed phase for large r(see Appendix 
B) 

a,o = £%)e-~^. (35) 

Q. 



q 



A similar behavior was found for a one-dimensional lattice in 



V. DISCUSSION 

A paired-fermion model for bosons with a hard-core interaction was studied. The effective 
Hamiltonian is given by Eq. ([1]). Both phases, the BEC and the MI phase, were found within 
the same mean-field approach to our model. 

In the BEC phase there is a long-range order in the phase fluctuations (i.e. phase coher- 
ence) and quasiparticles have a linear excitation spectrum for small momenta q. This reflects 
the existence of a Goldstone mode due to the breaking of a global U(l) gauge symmetry. 
Approaching the MI phase the system looses its phase coherence at the transition point and 
the excitation spectrum is characterised by the gap opening. 

In previous calculations, performed on the Bose-Hubbard model, each phase requires its 



own specific mean- field approach 20|, |21[ or a single one close to the phase boundary 8j. 
In the MI phase our expression for the excitation spectrum agrees with the branch of the 
excitation spectrum, which corresponds to creation of holes in the first Mott lobe of the 
Bose-Hubbard model, taking the limit of a large interaction. However, the second branch, 
which corresponds to formation of doubly occupied sites, does not exist in our model since 
we can create only holes in the singly occupied lattice to excite our system and not additional 
particles due to the hard-core condition. This is possible in the grand-canonical ensemble, 
where only the average number is fixed but the number of particles fluctuates. 
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Within a Bogoliubov approximation to the Bose-Hubbard model the quasiparticle spec- 
trum in the BEC phase was found as 



20 



e q = JJ 2 g 2 q + 2Un Jg 



where U is the interaction parameter and no is the condensate density. In contrast to this 
expression, we found for the spectrum the expression in Eq. (l2Tj) . These expressions do 
not agree in the limit U — > oo. Thus our hard-core Bose gas cannot be described within 
the Bogoliubov approximation to the Bose-Hubbard model by simply sending U to infinity. 
On the other hand, our results are in good agreement with a variational Schwinger-boson 
mean-field approach to the Bose-Hubbard model, which describe the phases near the phase 
transition, by sending U to infinity 

By changing the particle density from low to high values we pass from the weakly to 
strongly interacting regime. Within our approximation we were able to calculate an effect 
of quantum and thermal fluctuations on the condensate. The thermal fluctuations have 
strong effect: They destroy the phase boundaries in a direction of its vanishing. For fixed 
tunneling rate J and chemical potential u one can determine a critical temperature at which 
the condensate vanishes. The quantum fluctuations do not change phase boundaries. 



A. Critical Exponents 

Our expansion of the action A e g near the phase transition between the BEC and the Mott 
insulator Eq. (1231) is similar to the Ginzburg-Landau theory of (thermal) phase transitions if 
we use the relation 

-^f^— (36) 

and n c = J, i.e. the chemical potential u plays the role of the temperature T in our zero- 
temperature phase diagram. Within our approximation we find the mean-field expressions 
for three critical exponents: The critical exponent of the correlation length £ of density 
fluctuations is v — 1/2, of the compressibility \ it is 7 = 1 and of the order parameter 
it is f3 = 1/2. These critical exponents characterize the divergence of the corresponding 
quantities at the phase transition: 
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B. Static Structure Factor 



One possibility to detect the phase transition between the BEC and the MI in an optical 
lattice is to probe the excitation spectrum [l|, S]- I n the superfluid phase a broad continuum 
of excitations is observed but in the MI phase a more structured spectrum is measured, 
indicating the existance of the gap. The static structure factor vanishes in the Mott phase, 
as can be detected in experiments when studying phase transitions in an optical lattice. 

We also notice that the particle density cannot serve as an order parameter for the BEC- 
MI transition and the scaling law is not applicable for density-density correlations, which can 
be determined as a Fourier transform of the static structure factor and is given in Appendix 
B. However, we can assume the validity of the scaling law for correlations of the order 
parameter, namely 

(S(/>*Mr>). (37) 
Knowing the Green's function (see Appendix A) we obtain the following expression: 

(50 p (50_p') = Sp tP >G p ~ Spy—. (38) 

Then for large distances r and for J = fi (i.e. at the phase transition) we find the power law 

C , r ~ J d d r{5<Pp5^p)e lcir ~ J d d q^^^ r ~ JL , (39) 

which gives the anomalous exponent r\ = as in a mean field Landau- Ginzburg theory. 

VI. CONCLUSIONS 

We have used a paired-fermion model to describe strongly interacting bosons in an optical 
lattice with hard-core interaction. On the level of a mean-field theory we calculate the phase 
diagram, which includes the BEC and the MI. Including Gaussian fluctuations, we have 
found that the dispersion of quasiparticles is gapless in the BEC phase but has a gap in the 
MI phase: 

ef C = y/g^J* - ^ 2 ) + <7,V, 
e, MI = /x- J + Jg q . 

g q is the dispersion of the bosons on the lattice, defined in Eq. (jSJ). 

We have calculated the total density, the condensate density, and the static structure 
factor. We have shown that the quantum fluctuations as well as thermal fluctuations lead to 
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a depletion of the condensate, but the former do not change the critical points. The static 
structure factor contains information about the quasiparticle excitations and in the BEC 
phase for small q it is 

S{q) j 2 — coth— . 

It vanishes in the MI phase. For the critical behavior of the compressibiliy, the density- 
correlations and the order parameter at the phase transition between the BEC and the MI 
phase we found typical mean-field results. 



APPENDIX A: GREEN'S FUNCTION 



In this Appendix we write out the expression for the Green's function in both cases \(f>\ = 
and |0| 7^ 0. We denote p = {q, uo}. 



Case: \<j>\ = 



Deviation of the effective action due to fluctuations is 



5A cS = 50 p 5x P ) 



v- 1 - D{p) iD(p) \ I % 
iD{p) ± + D(p)) \5x P 



where 



D(p) = , , 1 , , v?= 1 



p Ji3-g q )' 
The determinant of the Green's function reads 



detGT' = %-Dip) 



(40) 



(41) 



Case: \<j>\ ^ 



Deviation of the effective action due to fluctuations is 



5A C{{ = ^( 5(p p ,5xp,S(p-p,Sx-p ) G 1 

p>0 



I Hp \ 

5(f)- p 
V $X-p J 



(42) 
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with the Green's function 

/ v- 1 - D(p) iD(p) 
iD(p) ±j + D(p) 



where 



G- l = 



—a 
ia 



la 
a 



\ 



\ 



—a 

ia 



ia 
a 



v- 1 - D(-p) iD(-p) 
iD(-p) ± + D(-p)J 



D(p) = 
D(-p) 



a = — 



1 \i 2 + J 2 + 2iiicu 

2 ' J(J 2 + uj 2 ) ' 

1 /i 2 + J 2 - %[IUJ 
~~ 2 ' J(J 2 + uj 2 ) ' 

1 |$| 2 /9 



2 J(J 2 + cj 2 )' 
The determinant of the Green's function is 



det G- 1 



1 



[2J 2 (3-^)] 2 (J 2 + ^ 



[u 2 + (f 2 -^) 9q + fiV q \. 



(43) 



(44) 



APPENDIX B: CORRELATIONS 



The decay of the density-density correlation function can be investigated as the inverse 
Fourier transform of the static structure factor 

(J 2 - /i 2 ) 



C r ,o 



J 2 



For large values of r the main contribution to the integral is for small values of q, namely 

1 



lim C, 



>0 ~ J dqdnq d e icF ~ 



— / 

-d+l I 



dqq sin(g) 



rd+l " 



(45) 



[1] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansen, and I. Bloch, Nature 415, 39 (2002). 
[2] D. M. Stamper-Kurn, A. P. Chikkatur, A. Gorlitz, S. Inouye, S. Gupta, D. E. Pritchard, and 

W. Ketterle, Phys. Rev. Lett. 83, 2876 (1999). 
[3] A. Gorlitz, A. P. Chikkatur, and W. Ketterle, Phys. Rev. A 63, 041601 (R) (2001). 
[4] D. M. Stamper-Kurn and W. Ketterle, Coherent Atomic Matter Waves, Les Houches 1999 

Summer School, Session LXXII p. 137 (2001). 
[5] L. Pitaevskii and S. Stringari, Bose-Einstein condensation (Claredon, Oxford, 2003). 

14 



[6] A. Y. Cherny and J. Brand, Phys. Rev. A 73, 023612 (2006). 
[7] C. Ates, C. Moseley, and K. Ziegler, Phys. Rev. A 71, 061601(R) (2005). 
[8] S. D. Huber, E. Altman, H. P. Biichler, and G. Blatter, Phys. Rev. B 75, 085106 (2007). 
[9] G. G. Batrouni, F. F. Assaad, R. T. Scalettar, and P. J. H. Denteneer, Phys. Rev. A 72, 
031601(R) (2005). 

[10] A. M. Rey, P. B. Blakie, G. Pupillo, C. J. Williams, and C. W. Clark, Phys. Rev. A 72, 023407 
(2005). 

[11] G. Pupillo, A. M. Rey, and G. G. Batrouni, Phys. Rev. A 74, 013601 (2006). 

[12] K. Ziegler, Laser Physics 15, 650 (2005). 

[13] C. Moseley and K. Ziegler, Laser Physics 15, 469 (2005). 

[14] J. Negele and H.Orland, Quantum Many-Particle Systems (Addison- Wesley, New York, 1988). 
[15] R. P. Feynman, Phys. Rev. 94, 262 (1954). 

[16] D. Jaksch, C. Bruder, J. I. Cirac, C. W. Gardiner, and P. Zoller, Phys. Rev. Lett. 81, 3108 
(1998). 

[17] G. G. Batrouni, V. Rousseau, R. T. Scalettar, M. Rigol, A. Muramatsu, P. J. H. Denteneer, 

and M. Troyer, Phys. Rev. Lett. 89, 117203 (2002). 
[18] M. Rigol, A. Muramatsu, G. G. Batrouni, and R. T. Scalettar, Phys. Rev. Lett. 91, 130403 

(2003) . 

[19] V. N. Popov, Functional Integrals in Quantum Field Theory and Statistical Physics (Reidel, 
Dordrecht, 1983). 

[20] D. van Oosten, P. van der Straten, and H. T. C. Stoof, Phys. Rev. A 63, 053601 (2001). 

[21] J.- J. Liang, J.-Q. Liang, and W.-M. Liu, Phys. Rev. A 68, 043605 (2003). 

[22] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslinger, Phys. Rev. Lett. 92, 130403 

(2004) . 



15 




FIG. 1: A projection of the phase diagram of the Bose-Hubbard model in the vicinity of the point, 
where the two Mott lobes meet. \i and J are in arbitrary energy units after the projection. 




FIG. 2: Phase diagram for ksT = 0. [i and J are in arbitrary energy units. 
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FIG. 3: Phase diagram for ksT = 0.05J. \i and J are in arbitrary energy units. 
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FIG. 4: Condensed density. The solid, dotted and dashed lines show the mean-field result at 
ksT = 0, the influence of the quantum fluctuations at ksT = to the mean-field result and the 
mean-field result at ksT = 0.2J, respectively. \i and J are in arbitrary energy units. 
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FIG. 5: Static structure factor in the BEC phase in the dilute regime, fi = —0.99 J. J is in arbitrary 
energy units. The unit of length is the lattice constant, so k is dimensionless. 
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